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Abstract 

The aim of this paper is to analyze different regulation mechanisms 
in spatial continuous stochastic development models. We describe the 
density behavior for models with global mortality and local establish- 
ment rates. We prove that the local self-regulation via a competition 
mechanism (density dependent mortality) may suppress a unbounded 
growth of the averaged density if the competition kernel is superstable. 
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1 Introduction 

We will discuss some classes on interacting particle systems (IPS) located in 
the Euclidean space M d . The phase space of such system is the configuration 
space T = T(]R d ) on IR d . By definition, each configuration 7 e T is a locally 
finite subset 7 C M. d . So, due to the standard terminology, we will deal 
with continuous systems. Random evolutions of IPS are given by Markov 
processes on T. Between all such processes, one may distinguish a subclass 
of so-called spatial birth-and-death Markov processes. In these processes 
points of a randomly evolving configuration appear and disappear due to a 
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Markov rule (see (3.1) below). Particular types of spatial birth-and-death 
processes are motivated by several applications. For example, Glauber type 
dynamics for classical continuous gases belongs to this class [T], [13]. Another 
very essential source of such processes is given by individual based models in 
spatial ecology or agent based models in socio-economic systems, see, e.g., [1] 
and the references therein. In any case, a concrete form of birth and death 
rates in the stochastic dynamics should reflect a microscopic structure of the 
system under consideration. 

To describe the problem we are going to analyze, let us start with the 
simplest case of a pure birth stochastic Markov process. In this process, new 
points appear in the configuration independently of existing points and loca- 
tions of these new points are uniformly distributed in the space. A possible 
interpretation of such random evolution is related to an independent creation 
of identical economic units in the space without any influence of their spatial 
locations. We will call this process the free stochastic development model. 
Another motivation to study such stochastic evolutions comes from appli- 
cations of the free development dynamics to generalized mutation-selection 
models in mathematical genetics. In these models a configuration describes 
locations of mutations inside of a genom and new mutations spontaneously 
appear and are equally distributed in the genom, see [22] , [16] , [8] . Obviously, 
this process is monotonically growing and it is easy to show that the density 
of particles in such a system will linearly grows in time. We would like to 
answer the following question: How may global regulations and local interac- 
tions change the asymptotic behavior of the system? More precisely, we will 
consider three particular cases of stochastic development models including: 

(i) A global regulation via a mortality rate that prescribes to particles 
(economic units) random life times (exponentially i.i.d. with a parameter 
m > 0). This case corresponds to the well-known Surgailis independent 
birth-and-death Markov processes on T, see, e.g., [23], [21], [H]. In the 
framework of mathematical physics, it is just Glauber dynamics for classical 
free gas. 

(ii) An establishment effect. In this case, the distribution of the position of 
a new particle depends on the local structure of the configuration. Newborn 
units will appear with small intensity in densely occupied regions. We will 
see that the establishment itself is not enough to prevent the growth of the 
density in the system. In fact, the establishment effects slower (logarithmic) 
growth contrary to the linear growth in the free model. 

(iii) A self-regulation via competition. The competition is reflected in a 
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density dependent mortality. The latter means that the mortality rate for 
each unit depends on the local structure of the configuration around this 
unit. The mortality rate enters into the model as a relative energy of a par- 
ticle inside the configuration corresponding to a competition potential. The 
described competition mechanism provides only a local regulation of dense 
regions inside a configuration. Nevertheless, Theorem 4.2 shows a global 
bound for the averaged density of the stochastic development process. Note 
that for the proof of this result we use an assumption of positive definite- 
ness (and, as a consequence, superstability) of the competition potential in 
the form stated in [19J. Therefore, the main result concerning the competi- 
tion case may be read as follows: a properly organized competition in the 
stochastic development systems produces a self-regulation for the density of 
units. 

Let us stress an essential point concerning the main aim of this paper. 
At present, we have quite restrictive conditions for the existence of general 
spatial birth-and-death processes, see e.g. J6]. In many applications we need 
a weaker information. Namely, we are interesting in the existence of Markov 
functions corresponding to given birth and death rates and a certain class of 
initial distributions on T. These Markov functions and their one-dimensional 
distributions are enough to describe the time evolution of initial states of sys- 
tems and to analyze asymptotic properties of stochastic dynamics (invariant 
states, ergodicity etc.). In the case of infinite particle systems, the concept 
of Markov functions is strictly weaker than that of Markov processes, and for 
particular models considered below there exist constructive methods which 
solve the existence problem, see [I], [6], [10] and [12]. But the main aim of our 
analysis is to obtain an a priori information about the time-space behavior of 
such important characteristics of these processes as correlations functions of 
their one- dimensional distributions which are probability measures on T. In 
particular, we are interested in the behavior of the particle density in course 
of the stochastic evolution. A constructive possibility to obtain some a priori 
bounds on characteristics of Markov dynamics may be also realized in other 
interesting IPS. As an example, we can mention the Dieckmann-Law model 
in spatial ecology where the existence problem remains open but conditions 
for explosion and non-explosion are stated in terms of the parameters of the 
model in [3]. Moreover, analogously to the well-known situation in the PDE 
theory, a priori bounds may play a crucial role in the study of the existence 
problem. 
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2 General facts and notations 

Let B(M. d ) be the family of all Borel sets in R d and Bb(WL d ) denotes the system 
of all bounded sets from B(R d ). 

The space of n-point configuration is 

rt } = rg d := { V C R d \ \v\=n}, nGN :=NU {0}, 

where \A\ denotes the cardinality of the set A. The space 

for A G Bb(M. d ) is defined analogously to the space . As a set, is 
equivalent to the symmetrization of 

= { (xi, . . . , x n ) e (R d ) n \ x k ^ Xl iik^l}, 

i.e. (R d ) n /S n , where S n is the permutation group over {l,...,n}. Hence 
one can introduce the corresponding topology and Borel <7-algebra, which 
we denote by 0(T^) and B(T^), respectively. Also one can define a mea- 
sure m*™) as an image of the product of Lebesgue measures dm(x) = dx on 
(R d ,B(R d )). 

The space of finite configurations 

To := U ^ 

neNo 

is equipped with the topology which has structure of disjoint union. There- 
fore, one can define the corresponding Borel a-algebra B(Tq). 

A set B G B(T ) is called bounded if there exists A G B b (R d ) and iVeN 
such that B C |J^=o ■ The Lebesgue — Poisson measure A z on T is defined 

as 

^-^ n\ 

?1=0 

Here z > is the so-called activity parameter. The restriction of A 2 to Ta 
will be also denoted by A 2 . We denote A := Ai. 
The configuration space 

r := {7 C R d I |7 n A| < 00, for all A G B b (M. d )} 
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is equipped with the vague topology. It is a Polish space (see, e.g., [9]). The 
corresponding Borel a-algebra B(T) is defined as the smallest a-algebra for 
which all mappings A^a : T — * N , N\(j) := (7 n A| are measurable, i.e., 

B(T) = a(N A \AeB b (R d )). 

One can also show that T is the projective limit of the spaces {T \} 
w.r.t. the projections p A '■ T — > Ta, Pa (7) := 7a, A G £>{,(!R d ). 

The Poisson measure n z on (T,B(T)) is given as the projective limit of 
the family of measures {ir^ }AeB b <u d )i where 7r^ is the measure on Ta defined 
by tt^ := e~ zm ^\ z . 

We will use the following classes of functions: L^Pq) is the set of all 
measurable functions on T which have a local support, i.e. G G L^Pq) 
if there exists A G £>&(IR d ) such that G [r \r A = 0; B^ S (T ) is the set of 
bounded measurable functions with bounded support, i.e. G |r Q \B= for 
some bounded B G B(T ). 

On T we consider the set of cylinder functions jFL°(r), i.e. the set of 
all measurable functions G on (T,B(T))) which are measurable w.r.t. B\(T) 
for some A G B b (R d ). These functions are characterized by the following 
relation: ^(7) = F |> A (7a)- 

The following mapping between functions on r , e.g. L^Fq), and func- 
tions on T, e.g. ^-"L°(r), plays the key role in our further considerations: 

KGW^^Gfr), (2.1) 

»?<E7 

where G G L^ro), see e.g. [3, [T71 HE]- The summation in the latter expres- 
sion is taken over all finite subconfigurations of 7, which is denoted by the 
symbol n d 7. The mapping K is linear, positivity preserving, and invertible, 
with 

K- X F{ V ) := ^(-l)l^lF(O, V G r . (2.2) 

Let A^f m (r) be the set of all probability measures /1 on (T, B(T)) which 
have finite local moments of all orders, i.e. J r \ , yx\ n p{d'y) < +00 for all 
A G Bb(R d ) and n G No- A measure p on (r , B(T )) is called locally finite 
iff p{A) < 00 for all bounded sets A from B(T ). The set of such measures is 
denoted by A^if(r ). 
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One can define a transform K* : A4f m (r) — > A4n(T ), which is dual to 
the i^-transform, i.e., for every /i G .Mf m (r), G G £>bs(ro) we have 

KG(i)n(d'y)= [ G(r,)(K*n){dr,). 
r Jr 

The measure p M := K*fi is called the correlation measure of /x. 

As shown in [7J for /x G Ai} m (T) and any G G L 1 (ro,p M ) the series (12. lft 
is /i-a.s. absolutely convergent. Furthermore, KG G L 1 (r,/x) and 



G(rj) p^drj) = /(JfG)( 7 )A*(*y)- ( 2 -3) 
r -/r 

A measure /x G A4f m (r) is called locally absolutely continuous w.r.t. n z iff 
/ia := ^°Pa 1 is absolutely continuous with respect to for all A G BA(K d ). 
In this case p M := K*fM is absolutely continuous w.r.t A. We define correlation 
functional k^ : T — > R + corresponding to the measure /x as the Radon— 
Nikodym derivative: 

The correlation functional k^ may be considered as the system of correlation 
functions corresponding to the restrictions k^ \ r (n). These functions are 
defined as following 

k (n) . ^dy „ R+ 

U n )(n r ) — I k fi ({x 1 ,...,x n }), if (xi, . . . , x n ) G (R d ) n 

(J, otherwise 



and they are well known in statistical physics, see e.g [20], [21]. In appli- 
cations a specially important role play correlation functions of the first and 
second orders: k^ l \x) and k^ 2 \x,y). These functions describe, respectively, 
the density of particles and pair correlations. 

A measure /x G A4} m (T) is called translation invariant if it is invariant 
with respect to shifts of configurations r 9 7 h- > {1 + a | a; 6 7} € T for any 
a G M d . The first-order correlation function of such measure doesn't depend 
on the space coordinate: k^(x) = fcjSr; and the second-order correlation 
function depends on difference of coordinates: kj?\x, y) = k^\x — y). 
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3 Stochastic development models 

Spatial birth-and-death processes describe dynamics of configurations in M. d 
when particles (agents, companies, economic units) disappear (die) from con- 
figurations and, on the other hand, some new particles appear (born) some- 
where in the space. The generator of spatial birth-and-death dynamics is 
heuristically given on measurable functions F : V — > R by 

(LF) (7) = d ^ T \ x ) \ x ) - 

+ / b(x, 1 )[F( 1 Ux)-F( 1 )]dx, (3.1) 

J«. d 

where d, b : M. d x T -> [0, cx)] are measurable rates of death and birth respec- 
tively. Of course, these rates should be finite for a.a. 7 G T with respect 
to a proper measure. Suppose that, additionally, b(-, r y) G Lj oc {R d ) for a.a. 
7 G T. Then this operator is well-defined at least on jFL°(r). Indeed, for 

f g ^L°(r) 

F( 7 \ x) - F( 7 ) = F(7 U x) - F( 7 ) =0, x G A c := R d \ A, 

and the both terms of (13.1 j) are finite. 

Note that if L is a generator of such process (even if we know that this 
process exists) then for the study of properties of the corresponding stochastic 
dynamics we need some information about the semigroup associated with L. 
This semigroup determines a solution to the Kolmogorov equation which has 
the following form: 

^ = LF t , F t U= F . (3.2) 

In various applications the evolution of the corresponding correlation func- 
tions (or measures) helps already to understand the behavior of the process. 
The evolution of the correlation functions of the process is related to the 
evolution of states of the system. The latter evolution is given as a solution 
to the dual Kolmogorov equation: 

— = LtH, [i t \ t= = fJ,o, (3.3) 

where L* is the adjoint operator to L on 7Wj m (r), provided, of course, that 
it exists. 
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Using explicit form of L we derive the evolution equation for quasi- 
observables (functions on To) corresponding to the Kolmogorov equation 
(13. 2p . It has the following form 

dG — 

-jf = LG t , G t \ t= =G . (3.4) 

Then in the way analogous to those in which the equation (13.31) was deter- 
mined for (I3.2p . we get an evolution equation for the correlation functions 
corresponding to the equation ( 13 .4ft : 

^ = L*h, k t \ t= =k . (3.5) 

The generator L* here is the dual to L w.r.t. the duality given by the following 
expression: 

((G, k)) = [ G-kdX. (3.6) 



Free development model A simplest model in considered framework is 
a model of free development when particles (identical economic units) are 
born independently without any influence of existing ones. An interpreta- 
tion is that a "decision" about appearing of a new company is produced 
outside of the system and it is not motivated by the situation inside of the 
system. Moreover, in this simplest model particles (units) will not die (will 
not become bankrupts). 

The formal pre-generator of the Markov dynamics that describes such 
model is the following: 

(L a F) (j) = a f [F(jUx)-F (j)} dx, 

where a > is the intensity rate of new units creation. 

It is easy to see that the operator L a is well defined, for example, on the 
set jFL°(r). The corresponding Markov process exists due to, e.g., [6]. Using 
results from [5j, we obtain 



(L a G) (rj) = a / G(f)Ui) dx. 



(3.7) 
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and 

(L:k)( V )=aJ2k(v\x). (3.8) 

x£rj 

Immediately from (I3.8P and (13.51) (see also [2]) we conclude that the den- 
sity of the free development model has the form 

k^\x) = k£\x) + at. 

Therefore, the density has linear growth in time. To prevent this growth we 
need to modify the generator introducing some regulation mechanisms in the 
model. 



Development model with global regulation Below we consider a model 
with a global regulation reflected in the death rate by an assumption about 
a finite life time for economic units. More precisely, we assume that each 
point of the configuration has exponentially distributed (with some positive 
parameter m) random life time and these random times are independent. 
Hence, a death (bankruptcy) appears as a random event equally distributed 
for all economic units independently of their space locations. 

A pre-generator describing such process has the following form: 

(L a> m F) ( 7 ) = m Y) [Hi \ x) - F( 7 )] + a I [F ( 7 U x) - F ( 7 )] dx. 

Note that the expression for L CTj m coincides with the one for the generator of 
so-called Surgailis process (see [23], [21], [H]). Again, using results from [5], 
we obtain 

(L CT , m G) (77) = -m\ri\G{rj) + a G(r]Ux)dx, (3.9) 
(L^ m k) ( V ) = -m\ V \k( V ) + (3.10) 

The considered stochastic dynamics has a unique invariant measure which 
is the Poisson measure on T with constant intensity — . 

Using (I3.10p (see also [2]) one can obtain a precise expression for the 
density of the process: 

kl 1 \y)=e- tm k^ (y) + £ (l - e+») . 
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Therefore, for an initial state with bounded density any positive global reg- 
ulation rate m gives time-space bounded density which converges uniformly 
in space to the limiting Poisson density. (For properties of higher order 
correlation functions see [2].) 



Establishment effects in the development model As we pointed out 
before, in the free development model an appearing of a new unit on the 
market and its location are independent of the presented configuration of the 
system. A reasonable generalization of this model is such that the newborn 
unit prefers to choose a location with smaller density of already existing 
units. The latter may be considered as higher probability to survive in less 
occupied regions. The corresponding term which decreases the birth rate of 
the generator in densely populated areas is called the establishment term. 
We consider the special case when this rate has exponential form, but our 
considerations may be extended to more general establishment rates. 
Let < G L 1 (R d ), <f>(-x) = (p(x), x e M d , and 

(VF)( 7 )= / exp{- S2(f>{x-y)}[F(rfUx)-F^)]dx. 

We suppose that there exist the dynamics of measures m and let k t will 
be corresponding correlation functions. Actually, the existence of a Markov 
process for considered case may be obtained from [6]. Using for any (p e 
C {M. d ) the equality 

9 C 9 



9t J Rd ^ V ^ dx= dt J ^' 7 ^ d ^ 1 ^ = ' L(t> ^' 7 ^ dfIt ^ 



r 



we obtain, by Jensen's inequality, 

= J exp|-J^0(x - y)}dfit (7) 




j/67 

exp 



- / (x - y) k { t l) (y) dy 



In the translation invariant case we obtain 

jk? > exp (- (0) k 



t / j 
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where (0) = J Rd (f)(x) dx. Hence, if g t is a positive solution of the equation 

= 

then /cq 1 ' ) > go implies fc$ > g t . 
One has 

e^dgt = dt 

J_ e {4>)gt - f + J^_ 
(0) " 

& = In «0>t + C), C>1 
(0) 

#o = ttt InC. 

w 

Putting for any feg the initial value go = fcg we obtain that 



A^^ln^t + expj^ (0) 



(0) 

Therefore, the establishment term cannot prevent unboundedness of den- 
sity. We may expect only essentially slower growth due to the establishment 
effect. 

Remark 3.1. Of course, if we consider two regulation mechanisms, namely, 
the global regulation and the establishment together, then the first- order corre- 
lation function will be also bounded (more precisely, all correlation functions 
will have sub-Poissonian bounds, cf. JM/). Moreover, in this case the operator 

(L G F) (7)=m^[F( 7 \i)-F( 7 )] 

£•£7 

+ a [F (j U x) - F (7)] expj- ^ <f>(x - y)}dx 

is the generator of so-called Glauber dynamics for a classical gas model (see, 
e.g., f77]/ . Jffl$ ). If (p has some additional properties such that there exists 
Gibbs measure with this potential, then such measure will be invariant (and 
even symmetrizing one) for the generator Lq. On the other hand, known 
properties of the corresponding Markov dynamics imply that corresponding 
correlation functions satisfied so-called generalized Ruelle bounds (see IT0{ 
for details). 
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4 Stochastic development models with com- 
petitions 

In the previous section we considered, in particular, global (outward) reg- 
ulation in the model. As wee see, such regulation may prevent unbounded 
(linear) growth (in time) of the density of our system. In this section we 
consider the case of a local regulation which appear due to the competition 
between elements (units) of the system. A pre-generator which describes 
such model has the following form: 

(La^F) (7) = J2 ( E a ( x ~ V)) [F(l\x)-F (7)] 



+ a [F (7 U x) - F (7)] dx. 

JR d 

Here < o 6 L 1 (R d ) is an even function s.t. 

(a) := a (x) dx > 0. 

J~R d 

The question about existence of a process with the generator L 0)CT we will 
not discuss in this paper. We just assume that there exist the dynamics of 
measures fit and let k t will be the corresponding correlation functional. 

Using results from [5] we obtain that 



(L a , a G) ( V ) = -2E a ( V )G( V )-J2(j2 a ( x - v)) G ^\*) 
+ a G(r]Ux)dx 

and 

(L* a a k) (77) = -2E a (t)) k (ry) - / V a (x - y) k (77 U x) dx (4.1) 



x) , 



xS7 
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where we used the following notations for the energy functional corresponding 
to the pair potential a(-): 

= ^2 a ( x ~ v)i v e r . 

{x,y}Cv 

It is easy to see that the Cauchy problems (13 .4p and (I3.5P for quasi- 
observables and correlation functions respectively have a form of hierarchical 
chains and, therefore, can not be solved explicitly. The latter is a common 
problem in the study of stochastic dynamics of IPS. In several particular mod- 
els such as Glauber type dynamics in continuum [T2], [TU] or some spatial 
ecological models [I] this difficulty may be overcame via a proper perturba- 
tion theory techniques. As a result, in the mentioned works we have existence 
results for corresponding evolutional equations together with certain a-priori 
bound for the solutions. Note that the perturbation techniques needs, in any 
case, a presence in the system a small parameter. In the considered model 
such parameter is clearly absent. Nevertheless, one can try to find estimate 
for the correlation functions. Actually, in the presented below approach we 
will use the explicit form of the Markov generator to obtain an a-priori bound 
on the density of the system. 

We will say that a sequence {A^, k £ N} of open bounded subsets of W d 
is of F-type if UfceN Afc = K , A& C Afc+i, k £ N and there exists F > such 
that for any h £ (0; 1) and for any k £ N 

\ A k\ 

where 



AJh) := < x : inf \x — y\ < h > . 

{ yeA k J 

A simple example of F-type sequence is the sequence of balls A& = B(0, k) 
with center at origin and radius k £ N. Indeed, for any h < 1 

For any A £ i3 c (lR d ) we will call the average density of the our system the 
following object 



|a[ J ^' \ x )dx, 
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where k^ is the first-order correlation functions (density) at moment t > 0. 

Lemma 4.1. Suppose that the function a is continuous and positive definite 
and the sequence {A k , k G N} is F-type. Then there exists c > such that 
for any open A G {A^, fcGN} 



12 

|aT 



2E a (r 1 )>c^- ) r]eT A . 



Proof. In [19], it was shown that for any continuous positive definite function 
a the energy E a is superstable, namely, for any open A C M d and for any 
rj := {xi}™ =1 C A the following inequality holds 

2E M >" 2 l{a) - Hh)]2 



where 



<\x\ 

Therefore, for any A e {A^, fcGN} 



|A| [(a) + 5 (h) + a (A, h) (a)]' 



5(h) = 2 i a (x) dx > 0. 

J\x\>h 



n 2 [(a) -5(h)] 



n 2 



2Ea{v) - \A\ [5 (h) + (F + 1) (a)} : |A| C ' 
Let h G (0; 1) be such that 

(a) — 5 (h) = / a(x) dx — a (x) dx ^ 

J\x\<h J\x\>h 

(we may always choose such h since the first integral is an increasing function 
of h and the second one is a decreasing function). Stress that c > and 
doesn't depend on A. □ 

Theorem 4.2. Suppose that the function a is continuous and positive definite 



and the sequence {A^, k G N} is F-type; let c be as in Lemma \4-1\ Suppose 
also that the: 
t > 0, k G N 



also that there exists D > \j — such that p$ k < D, k G N. Then for any 



pt k < D. 
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Proof. Note that for F (7) = (<p, 7), 7 G T, <p G Co(Mr) we have 

{L a ,*F) (7) = - E a (x - y> (x) + <r / ^ (x) da:. 

Let (a:) = 1 A (a:), A G {A fe , fc G N}. Then F( 7 ) = | 7a | and 



(l 0iCT f) (7) = - e E a (*-tf) I ^W+H 

= "E f E a(x-y)\+a\A\ 

<-E ( E a ^-^)) +^i a i 

= -2£ a ( 7A ) + a |A| < -t^t | 7 a| 2 + o- 



|A| 



Let us set 



nf := jf | 7 a| dtni-y) = J (1a, 7) <Wt) 



l A (x)£; t (1) (x)dx = / ^ iJ (x)dx= |A|pf. 
' A 



Then using Holder inequality 



= jf L a a | 7A | (7) < I (a I A| - JL | 7A | 2 ) (7) 

= * |A| - ^ jf | 7A | 2 <fo (7) < a I A| - A Qf | 7A ] dn t (7)) ' 



-IAI-t^K) 2 . 



As a result, 



A \2 



S rf<-c(rf) 



16 



D. Finkelshtein, Yu. Kondratiev 



Therefore, if we consider the positive solutions of the Cauchy problem 

j t 9(t) = a-cgHt) (42) 
9(0) =9o 

with proper g > and if p§ < go then pf < g (t), t > 0. Solving (14. 2f) we 
obtain 

ln |^ + ^| -lnC = 2v /^ C>0; 

Then 

«(°» = V!( 1 + c^t)' CeR ' 

Let go = D > Pq . Then since -D > we have 

C= = + 1 > 1 

and 

Ce 2 ^* - 1 > C - 1 = = > 0. 

As a result, 



for any t > and for any A G {A&, k G N}. The statement is proved. □ 

Corollary 4.3. Under conditions of Theorem \4-S\ in the translation invariant 
case we have that k^ < D implies k[ < D. 

At the end we consider a simple estimate for the second-order correlation 
function. Let us suppose that 

a(u)>0, ueR d . 



Regulation mechanisms in spatial stochastic development models 



17 



Then in the translation invariant case the following estimate holds 

kf ] (u) < e- 2a ^ l k^ (u) + 2a f e^^k^dr 

Jo 

< e- 2a{u)t k {2) (u) + 2aD [ e^^^dr 

Jo 

= e- 2a ^k^ (u) + (1 - e- 2a(u)i ) 

We have two possible estimates 

kf\x -y)< e-^-yH^ (x - y) + ° D (4.3) 

a(x-y) 

and 

kf ] {x-y)< e- 2a{x - y)t k {2) (x-y) + CaDt. (4.4) 
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